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'■^ ' Abstract. The semi-classical stability of several AdS NUT instantons is studied. 

i . Throughout, the notion of stability is that of stability at the one-loop level of Euclidean 



Quantum Gravity. Instabilities manifest themselves as negative eigenmodes of a modified 
■ Lichnerowicz Laplacian acting on the transverse traceless perturbations. An instability is 

found for one branch of the AdS-Taub-Bolt family of metrics and it is argued that the other 
' branch is stable. It is also argued that the AdS-Taub-NUT family of metrics are stable. A 

j^. . component of the continuous spectrum of the modified Lichnerowicz operator on all three 

CNj | families of metrics is found. 
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O ■ 1. Introduction 

Motivated by the AdS/CFT conjecture Q, there has recently been much interest in the 
two parameter AdS-Taub-NUT familyf of Riemannian biaxial Bianchi-IX metrics satisfying 
the Einstein equations with negative cosmological constant and non-trivial NUT charge 

El SI 13- 

Upon the imposition of a regularity condition such metrics divide into two one parameter 
classes. The first class (AdS-Taub-Nut) have self-dual Weyl tensor and contain a nut. 
The second class of solutions (AdS-Taub-Bolt) contain a bolt, this class splits further into 
two branches. This is analogous to the case of AdS-Schwarzschild solutions at a given 
temperature, where the role of AdS is played by the AdS-Taub-Nut. The AdS-Schwarzshild 
solution with a smaller mass (and smaller horizon) is unstable and it's action is greater than 
that of both AdS and the other AdS-Schwarzshild MM- We find a similar situation for AdS- 
Taub-Bolt. 

In 0, Hartnoll and Kumar conjectured, based on the Klebanov-Polyakov version of 
the AdS 4 /CFT 3 correspondence, that the global minimiser of the action for the AdS-Taub- 
NUT class of metrics should be stable. In ifTOl . the stability of some of these spaces against 
scalar perturbations and brane nucleation was discussed. By considering a one-loop correction 
to the bulk gravitational partition function, we shall investigate the semi-classical linear 
stability of such spacetimes, using techniques developed by Hu [11] and applied to the case 

| Following |2), we refer to the full two parameter family as AdS-Taub-NUT, reserving AdS-Taub-Nut for the 
regular solution containing a nut 
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of Euclidean Taub-Bolt by Young lTT2"ll. These techniques have also been applied to the case 
of Lorentzian Taub-NUT recently by Holzegel 1131 . The criterion for instability is that the 
modified Lichnerowicz operator 

A L h ab = -V e V e /i ab - 2R acbd h cd (1) 

acting on transverse, trace-free symmetric tensors should have no negative eigenmodes. We 
shall investigate the spectrum of this operator on metric perturbations of the AdS-Taub- 
instantons which preserve the SU (2) symmetry. We find a negative mode for one of the 
instantons indicating instability for all values of the cosmological constant. We also find part 
of the positive spectrum of the modified Lichnerowicz operator for all instantons in the class 
under consideration. 

In Section 2, we give a brief overview of the AdS-Taub-Nut and AdS-Taub-Bolt spaces, 
stating some results which we will later require. In Section 3, we introduce the method used 
to examine the stability of the spaces and we present an instability for one branch of the AdS- 
Taub-Bolt spaces and argue that the other spaces are linearly stable. In section 4 we present 
an alternative viewpoint, confirming the claims of section 3. In Section 5 we investigate the 
continuous spectrum of the operator A^. 



2. NUTs and Bolts in AdS space 

The metric for both AdS-Taub-Nut and AdS -Taub-Bolt in four dimensions can be put into the 
form: 

ds " = -TT^ dr2 + ^ 2 A(r)a 2 + B(r)(a 2 + a 2 2 ), (2) 
Ayr) 

where cr^ are the left-invariant £77(2) one-forms. The functions A(r) and B{r) are given by 

r 2 + N 2 — 2mr + r 2 (r 4 - QN 2 r 2 - 3iV 4 ) 



A(r) 



r 



2 



N 2 



B(r) = r z - N\ (3) 

where N is the NUT charge and m is the mass. This metric is Einstein, with = —^g^ u . 
The fixed point set of the U(l) action is given by the solutions, r = r + of A(r) = 0. If 
B{r + ) = 0, the fixed point set is of zero dimension and is known as a nut. If B(r + ) ^ 0, 
the fixed point set is two dimensional and is known as a bolt. As r approaches r + , we will in 
general have a conical singularity unless A(r) satisfies the regularity condition 

l^(r+)l = 2F W 
This amounts to a relation between m and N, which can be solved and we find that the number 
of solutions depends on the value of jr. 

For all values of N and £ there is a solution, called the AdS-Taub-Nut solution where 

N(£ 2 - 4A 2 ) 

r n = N, m n = V '-. (5) 
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Clearly we have B(N) = 0, so we have a regular nut at r = r n . We can also find two bolt 
solutions, which I shall call the AdS -Taub-Bolt ± solutions. They are given by: 



£ 2 ± v^ 4 - 48iV 2 £ 2 + 144A 4 

n± - ' 

r 4 ± + (£ 2 - QN 2 )r 2 b± + N 2 (£ 2 - 3A 2 ) 

mb± = 2^ • (6) 

Clearly these equations only make sense if the quantity under the square root in the first 
equation is positive. This restricts £/N to the range 



— > 2^3(2 + ^) ^6.69. (7) 

Thus we find that requiring that the metric © be regular restricts the freedom quite 
considerably. For any given value of -4, there are either 1 or 3 regular metrics in this family, 
excluding the critical case where r b+ = r b _. It is possible to calculate the action for these 
spaces, but we shall postpone this until section 4. 

It is useful to consider the limit i — > oo, as we expect to recover the Euclidean Taub-NUT 
and Taub-Bolt solutions. For the AdS-Taub-Bolt + case, this limit does not exist as r b+ — > oo 
as £ gets large. However, for the other two cases we can take this limit and we find for the nut 
case: 

r - N 

m n -> N, A{r) -> — — , (8) 

r + JM 

which gives the well known form of the metric for the self-dual Taub-NUT instanton. For the 
AdS-Taub-Bohr case, we find: 

5 r 2 - %Nr + N 2 
r b _->2iV, m b _^-N, A{r) , (9) 

which gives the metric for Euclidean Taub-Bolt. This will provide a useful check on our 
stability results as for vanishing A, Taub-NUT has a self-dual Riemann tensor, and hence is 
linearly stable lfT4l . In |Q2J an unstable mode for the Taub-Bolt instanton was found. 

3. Stability 

3.1. A criterion for instability 

Instabilities of a physical system may make themselves known as an imaginary part of the 
partition function for the system lfT5l fT6l . The partition function for Euclidean quantum 
gravity is given by: 

"d[g]e- s ^, (10) 

where the integral is taken over all Riemannian metrics subject to appropriate boundary 
behaviour and periodic in Euclidean time. The Euclidean action is given by: 

S[g]- 



[ dVV5(i2-2A)--L ! dS^/hK. (11) 
167TG J M 8ttG J 9m 
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Unfortunately this is not positive definite. Under a conformal transformation of the metric, 
g' = Vt 2 g, where for simplicity we assume Q — 1 on a neigbourhood of dAi the scalar 
curvature transforms like: 

R! = n~ 2 R - 6fT 3 V a V a ll (12) 

Thus we find 

J [s1 = ~TiT~n I dV ^ {&R + QQ a ^ a ) - -L / dSVhK. (13) 
lb7rG J M 87TG J dM 

So by taking f2 to vary quickly, we can make I[g'] as negative as we choose. This problem can 
be circumvented by an appropriate choice of contour, but only in the semi-classical (one-loop) 
approximation ifTTIl . For a semi-classical approximation, we expand the integral (flOb about 
the critical points of S[g], where we expect the dominant contributions to the integral. These 
satisfy 

6S 



0^R ab = Ag 



ab ■ 



vg a b 

i.e., the classical Riemannian vacuum Einstein equations. Thus in order to find the semi- 
classical approximation to Z, one expands in small perturbations about the classical solutions: 

Slab = gab + h a b, 

S[g] = S [g] + S 2 [h]+O(h 3 ), (14) 

where S 2 [h] is second order in the small perturbation /i afe . We truncate the series for I[g] and 
integrate over h ab to get 

Z = Z lloop = Ne- 30 ^ j d[h}e~ s ^ h \ (15) 

where the integral is to be taken over physically distinct perturbations h ab . The action S 2 [h] is 
invariant under gauge transformations which correspond to infinitesimal diffeomorphisms: 

K b = h ab + V a V b + V b V a . 

In order to deal with this, we follow ifTHI and use the Fadeev-Popov gauge fixing technique. 
The gauge fixing condition will be: 

V tt f h ab - \g ab h] = W\ (16) 



P 

where W is an arbitrary vector, (3 is an arbitrary constant and h = h a a is the trace of the 
perturbation. The standard Fadeev-Popov method allows us to re-write the functional integral 
(1751) as an integral over all fields h ab , but with an altered integrand: 

Ziioop = Ne~ So[9] J d[/i](det C)e- Ja « ff - [h] . (17) 

We now decompose the metric perturbation into a transverse, traceless part 4> ab , a longitudinal 
traceless part generated by the vector h and the trace h. Further, a Hodge-de Rham 
decomposition can be performed on h so that 

n a = Va + V a x, with V a r/ a = 0. 
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We can then write the effective action as: 

1 



'2eff. 



dVC 2 , (18) 



16vr 
with 

£ 2 = -l^Gabdt * + \hFh - 2 1 D a \x)D ab {Fx) - l^ abe (Ve)^a b f((C 2 y d Vd ) (19) 
where 7 is a constant related to (3 by (3 = 47/(1 + 7) and we have defined operators 

F = l^v e v e - h, 

I67 4 ' 

D ab = 2V a V fe - -g ab V e V e , 

T> abe {y e ) = \{v a v b + v b v a ), 

(C 2 ) ab V b = - (V e V e + A) V a , 

G abcd ^ cd = - V e V e afe - 1R acbd ^ c \ (20) 

where V is a divergence free vector. It is possible to write the (det C) factor in (fT7t as a 
functional integral over anti-commuting vectors. Applying a Hodge-de Rham decomposition 
one finds 

(det C) ~ (det C 2 ) (det F). 
Finally, evaluating the Gaussian integrals in ( fT7l we find 

Ziioo P ~(detG)^(detC 2 )i (21) 

Any zero modes of the determinants in (EH should be projected out. The determinants can be 
regularised by a ^-function regularisation lfl~8l . C 2 is a positive semi-definite operator on the 
space of divergence free vector fields, with zero modes corresponding to Killing vectors. G 
however is not positive definite and can in fact have a negative eigenmode. Such a negative 
eigenmode would introduce an imaginary part to the partition function and thus herald an 
instability. For example the flat Schwarzschild solution has a negative eigenmode. It is to be 
expected that the partition function be pathological in this case since the canonical ensemble 
for black holes breaks down due to the fact that they have negative specific heat. 

Thus in our search for instabilities we can restrict to perturbations h ab which are 
transverse and tracefree. Our criterion for instabilities is that there exist negative eigenvalue 
solutions to the eigenvalue equation: 

(A L h) ab - 2Ah ab = -V e V e h ab - 2R acbd h cd = Xh ab , (22) 

where we have re-expressed G in terms of the Lichnerowicz Laplacian A L lfT9ll . This equation 
is consistent with the Transverse Traceless condition. 

A negative eigenmode of d22l corresponds to a direction within the space of gauge fixed 
metric perturbations along which the classical solution g ab is a local maximum. The criterion 
that 

At = At - 2A > 
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for stability may be arrived at from a more geometric perspective (see for example Besse lEUll ) 
which avoids the Fadeev-Popov procedure used above. There are some cases in which it can 
be shown that Al is positive definite. These include the case when g a t has self-dual Riemann 
tensor (sometimes called the half-flat condition) ifTHl and also the case where g a b is Einstein- 
Kahler ll2T1 . In both these cases, the existence of a covariantly constant spinor is used to relate 
the eigenmodes of Ax, to those of scalar Laplacians acting on charged fields. 



3.2. Hu 's Technique 

Hu ifTTIl proposed a method for separation of variables of tensor equations in a homogeneous 
space where the metric can be written in the form 

ds 2 = dr 2 +^f~ j a i a\ (23) 

where i,j = 1,2,3 and a 1 = Ui are again the left invariant one-forms on 577(2). This 
was generalised by Young [12J to the case where goo is permitted to depend on r. The 
metric perturbation is expanded in terms of the Wigner functions, V KM J (9, (ft, ip) which 
are the analogue on S 3 of the spherical harmonics Yi m {9, <p). We concentrate on diagonal§ 
metric perturbations with the lowest "angular momentum", J = M = K = 0. These are 
the perturbations which preserve SU (2) symmetry. The resulting equations, along with the 
transverse and traceless conditions for the general metric form © are given in the Appendix. 

We first consider the 11 and 22 equations. Taking the difference of (IA.1I) and (IA.2I) we 
find a second order differential equation for Y = h u — h 22 . We can put this into Schrodinger 
form: 

- + V ( r *)* = X *i 0<n< oo, (24) 

via the substitutions Y — fx and dr = gdr*, for some suitable / and g. When we do this, we 
find that V(r(r*)) is a positive function for r > r + , so that no normalizable solutions to d24l 
exist with A < 0. This means that in the search for a negative eigenmode of d22l we may set 

hu = h 22 . 

We now consider the other diagonal mode. Setting hu = h 22 , we can use the constraint 
equations (IA.5I) . (IA.6I) with the 00 equation (IA.3I) to decouple a second order differential 
equation in h 00 from the others. If we can solve this, we can find h n , h 22 and /i 33 from 
the constraint equations. Since d22l is consistent with the transverse traceless condition, we 
know that (IA.5I) . (IA.6I) and (IA.3I) imply (IA.1I) and (IA.4I) . This can be checked explicitly. The 
equation we decouple may be written: 

a(r)/i ' + b{r)h' m + c(r)h 00 = 0, (25) 

where differentiation with respect to r is denoted by a prime and the coefficients are given by: 

a(r) = A, 

AB' B(A' 2 -2AA") A f 2AB' 2 — BA'B' — 2ABB" 
b(r) - 3A + _ + _ _ _ - - ^ _ _ _ 

§ The equations for off-diagonal perturbations can be put into Schrodinger form with strictly positive potential 
(see i24\ ). for all three metrics considered here and so cannot give rise to negative eigenmodes. 
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, . , A 12 A'B' AB' 2 „„ 3 fBA' + AB' 
c ( r ) = A + 77T + ^ 55-+^ 



2A S 5 2 2A \BA' -AB' 

BA'B' + 25A5" - 2AB' 2 (IB A' + AS' 
+ B 2 V BA ' ~ AB ' 

the rest of this section will consist of an analysis of this equation. 



3.3. AdS-Taub-Nut 

In the case where A(r) and B{r) are given by © with © we can cast equation d25t into 
Schrodinger form d24l) where V is once again found to be positive on the range r > r+. Thus 
for AdS-Taub-Nut we have checked all the possible TT perturbations with J = K = M = 
and found that none give negative eigenvalue solutions to (|22b . thus there is no linear 
instability due to perturbations of this type. This is in agreement with the result that the 
£ — ► oo limit yields a (linearly) stable metric. 

One would normally expect the eigenvalues of a Laplacian operator to be bounded below 
by the most symmetric modes. In our case that would be the SU{2) invariant modes. In 
the case of the scalar Laplacian acting on a space with metric given by © this can be seen 
explicitly: 

_ V 2 0= - V 2 ^ K {r)V K {6^^)) 

= (- V 2 ^(r)) V K (6, <j>^)-2 (dU K (r)) (d d V K (9, 0, if;)) 

= (- V 2 ^(r)) V K {9, 0, V) + <t> K {r) (-V 2 V K (6, 0, . (26) 

It can be shown that (— V 2 V K ) is positive for this metric form. This separation into a radial 
part and an angular part does not occur for the the Laplacian acting on symmetric tensors as 
there are extra terms introduced by the connection. It has not been possible to show explicitly 
that the SU(2) invariant perturbations give a lower bound for A L , but one would still expect 
this to be the case. Increasing the spin of the perturbation modes introduces a more rapidly 
varying angular dependence, which one would expect to increase the eigenvalues of —V 2 and 
hence of A^. 

We therefore conjecture that A L is positive acting on all metric perturbations and thus 
that AdS-Taub-Nut is stable. This may be related to the fact that the Weyl tensor is self dual 
for the AdS-Taub-Nut space. In the Euclidean case, the full Riemann tensor is self-dual and 
this is known to imply stability iflH . 



3.4. AdS-Taub-Bolr 

In the case of AdS-Taub-Bohr, when one changes (l25t into Schrodinger form, one encounters 
a divergent potential. In order to proceed, we must use numerical integration techniques. If 
we make the substitution r = xN, we find that the equations depend on only through the 
ratio jr. We can therefore scale N out of the problem and set N = 1 without loss of generality. 
We must now consider the singularity structure of the differential equation (l25l) . In the region 
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of interest to us, r > r b _, we find that the equation has 3 regular singular points, at r b _, r int 
and infinity, where r b _ < r int < oo. Other singular points outside of the range of interest 
preclude the possibility of an analytic solution. 

We consider first the regular singular point at r = r b _. Near to this point, the differential 
equation has the asymptotic form: 

K + — — h' 00 + —^h m = 0. (27) 

r — r b „ (r — r b -J^ 

Substituting h 00 ~ (r — r b _) Q , we find the indicial equation 

a 2 + 3a + 2 = 0, 

with solutions a = — 1 and a = —2. Thus there are two independent solutions of (l25t that 
behave like (r — r^)" 1 and (r — r b _)~ 2 as r — > r b _. We impose a normalisation condition 
on the perturbation modes given by 

dVh ab h ab < oo. (28) 

Using the constraint equations, we can show that the h 00 contribution dominates the integrand 
in a region just outside the bolt. Using the fact that dV = ^/gdipd9d(j)dr, where xj), 9 and are 
the usual Euler angles on SU(2) and the form of the metric © with the constraint equations 
dA.5UA.6l) we can show that this integral coverges at r = r b _ if and only if h 00 ~ (r — r h _Y 
where j3 > —3/2. Thus only one of the two independent solutions at r = r b _ represents a 
normalisable perturbation - we must pick the solution that behaves like h 00 ~ (r — r b _) _1 
near r = r b _. This provides us with a stepping off condition for our numerical integration. 
We can perform a Frobenius expansion of the solution about the regular singular point and 
then use this to calculate h 00 and h' Q0 at r = r b _ + e where e is some small number, taken to 
be 10~ 5 . This deals with the first regular singular point. 

It can be shown, by a similar method to that used above that h 00 and h' Q0 are bounded as 
r r int.> so no special considerations are required for the numerical integration through this 
point. 

Finally we need to take account of the regular singular point at infinity. In the limit that 
r — ► oo, the differential equation has the asymptotic form 

Ko + -Ko + ^±r^oo = 0. (29) 
Substituting h 00 ~ r Q we find the indicial equation: 

a 2 + 11a + (28 + £ 2 X) = 0, 

with solutions 

-ll±V9-4l 2 A 

a± = . (30) 



Using the constraint equations, we find that the integrand of (1281) is dominated by the 
contribution from h 33 ~ r 6 h 00 . This translates to a requirement that h 00 ~ r 13 with j3 < — y in 
order that (l28l is satisfied. So we see that only one of the two independent solutions at infinity 
will give a normalisable perturbation. We will need to match the normalisable solution at 
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Figure 1. A plot showing the negative eigenvalue, A, against I for AdS-Taub-Bolt Space 



infinity to the normalisable solution at r b+ and this matching, if it is possible, will determine 
the value of A. 

In order to ascertain whether there exist negative solutions for A, we perform a numerical 
integration starting at r b _ +e, with initial conditions determined from the Frobenius expansion 
about rt,_. This determines the solution, up to a arbitrary multiplicative constant as we can 
determine h 00 and all its derivatives here. Having chosen the normalisable solution at the bolt, 
we need to check that the solution is normalisable at infinity. The only parameter we still have 
available is A. As r — ► oo the solution must have the asymptotic form found above, but will 
in general be a linear combination of the two independent asymptotic forms. 



where a + > —11/2 and a_ < —11/2 are given in equation d30b and k(X), /(A) are constants, 
which we assume to depend continuously on A. The normalisation condition requires that 
A = An, where k(X ) = 0. Since we are only performing a numerical integration, it is not 
possible to find the form of k(X) explicitly, but we can find an interval within which A must 
lie. 

If we consider f(r) = r ll l 2 h m {r), men f( r ) w iU genetically tend to either positive or 
negative infinity, depending on the sign of k(X). If we can find numbers Ai and A2 such that 
fc(Ai) has opposite sign to k(X 2 ), then we can deduce the existence of a root of k(X) in the 
interval (Ai, A2), by the Intermediate Value Theorem. This procedure was implemented using 
the computer package Mathematica, and negative eigenvalues have been found for values of 



£ on the entire range 2y 3(2 + y/3) < £ < 00. A plot of these negative eigenvalues is shown 
in Figure 1, together with a best fit curve. 

We find that as £ — > 00, we have that A — > —0.20 approximately. This is consistent 
with the findings of Roberta Young 1121 . who found that Euclidean Taub-Bolt had a negative 



h 00 (r) ~ k(X)r a + + /(A)r' 



(31) 
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eigenmode in this sector with A ~ —0.20. We also find that as t approaches the critical value 
where the AdS-Taub-Bolt and AdS-Taub-Bolt + spaces are the same, the negative eigenvalue 
tends to zero from below. 

3.5. AdS-Taub-Bolt + 

The final case, where the bolt is located at r = r b+ is similar in many ways to the previous 
case. We find the same singularity structure and asymptotic forms for the perturbation. We can 
proceed in exactly the same way as above, but we find no negative eigenvalues. We therefore 
conjecture that this spacetime is stable, as it appears to have no negative eigenmodes in the 
lowest "angular momentum" state. 




Figure 2. A plot showing the actions S — Snut against i for N = 1 



In Figure 2 we show a plot of the action of both AdS-Taub-Bolt spaces against i, with 
the action for AdS-Taub-Nut subtracted. This will be calculated in Section 4. The dotted line 
represents the unstable AdS-Taub-Bolt" space, the horizontal axis the Taub-NUT space and 
the other curve AdS-Taub-Bolt + . The thicker line denotes the global minimizer of the action. 

We see here two interesting features. Firstly we note that at £ = 2y / 3(2 + y/3) we have a 
bifurcation, with the two Bolt families appearing, with one stable and the other unstable. We 
also note that at £ = 2 a/7 + 2-\/T0 there is a phase transition from AdS-Taub-Nut as the global 
minimiser of the action to AdS-Taub-Bolt + as the global minimiser. This is the NUT charged 
version of the Hawking-Page transition for AdS Black Holes 0. 
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4. Another Perspective on Stability 

We will now present an alternative argument, which reproduces all the features of the plot in 
Figure 2. 

In constructing the nut charged AdS instantons considered above, the usual technique 
(outlined in Section 2) is to consider a metric of the form © which satisfies the Einstein 
equations = —j^g^u- Then one imposes the condition of regularity to reduce the number 
of solutions to those given above. Here we shall consider a family of metrics which are regular 
everywhere, with a bolt (or possibly nut) located at r b , which is treated as a free parameter. 
The AdS-Taub-Nut and Taub-Bolt(±)-AdS instantons correspond to particular choices of r&. 
We can then calculate the action for this family of spacetimes as a function of r b and we 
expect it will be locally extremized at the already known values of r b which give solutions 
of the Einstein equations. This corresponds to taking a slice parameterised by r& through the 
space of metrics containing a bolt or a nut. 



ds 2 = -r^dr 2 + AN 2 A(r)al + B(r)(af + a 2 ), (32) 



with 



Our metric ansatz will be: 
1 

A{r 

(r — r b )(ar + j3 + £~ 2 (r 3 + r 2 r b )) 



A(r) 



r 2 _ ]\[2 

B(r) = r 2 - N 2 . 

We then impose regularity at r = r b and require the mass to be m, by considering the 
asymptotics. This gives us the metric functions: 

(r — r b ) (4N (r — r b ) r b (r + r b ) 2 — I 2 (4mA (r b — r) + (r + rb) (n 2 — r 2 ))) 
A{j > = A£ 2 N (r 2 - A^ 2 ) r b ' 

B(r) = r 2 - N 2 . (33) 

We are now ready to calculate the action. The Euclidean action is given by: 

S = S hulk + Ss.rf. = / dV^ ("i? + A ) _ J_ f dSVhK, (34) 

16vrG J m V £ J 87rG JdM 

which is the standard Einstein-Hilbert action together with the Gibbons-Hawking-York 

boundary term. It is well known that this integral does not converge and so some means 

of regularising the action is required. The method used was that of counterterm subtraction 

proposed by Emparan, Johnson and Myers ll22ll . This was also found independently by Mann 

ll2"3l . The calculation proceeds much as their calculation for the action of Taub-NUT. We find 

that the action of the metric defined by d32l and d33t is: 

S=^-(r 3 b -3N 2 r b )+N 2 + 4Nr b -r 2 b . (35) 

We note that this does not depend on m, which is also a free parameter of this family of 
metrics. The variation of this reduced action alone will not determine m, for that we require 
the full Einstein equations, which force m to take the NUT or Bolt values at the critical points. 
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Setting r b = xA^defining ^ — £/N, and using the additive freedom to set S 
Taub-Nut we find 

-J- = - 3x + 2) - 3 + 4x - x 2 . 

N z fx 1 



S 




Figure 3. A plot showing S/N 2 against x for several values of fx = £/N 



So we see that as fx is varied, the graph of S/N 2 will change. A plot of S/N 2 against 
x for various values of fx is shown in Figure 3. Also included is a dotted curve showing the 
locus of the extremal points of the function as fx varies. We see that the point x — 1, which 
corresponds to AdS-Taub-Nut is always a minimum since we must exclude the region x < 1 
from consideration as these metrics are not regular. For small fx this is the only minimum. 
As we increase fx, we find that a pair of extrema, one maximum and one minimum appear at 
fx = fxi, where 

l x 1 = 2^3(2 + (36) 

which is precisely the value we would expect from the standard constructions for AdS-Taub- 
Nut and AdS-Taub-Bolt. If we think of this function S(x)/N 2 as defining the dynamics of 
some system by 

■± = -S'(x)/N 2 , 

then this bifurcation is of the form known as a saddle-node bifurcation, since it produces two 
new fixed points, one stable (a node) and one unstable (a saddle). Of course we have no 
formal dynamics on the space of metrics, but we still expect a local maximum to correspond 
to an unstable metric as it corresponds to a direction contributing an imaginary part to the 
partition function. Thus we shall refer to this bifurcation as a saddle-node bifurcation. We 
also find that the maximum occurs at precisely x = r h _/N and the minimum at x = r h+ /N, 
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confirming our previous result that we expect the AdS-Taub-Bolt instanton to be unstable 
while the AdS-Taub-Bolt + instanton is stable. 

We also note a global bifurcation when fi = ^ 2 where 



when the global minimum moves from x = 1 to x — r^ + /N, which corresponds to AdS- 
Taub-Bolt + . This again accords with what we expect from previous analysis and corresponds 
to the Hawking-Page type phase transition. 

These observations accord nicely with the general argument given in section 6 of ll24l 
for the production of negative modes of the Lichnerowicz operator at bifurcation points in 
parameter space, in particular, from figure 1 we find that the negative eigenmode tends to zero 
as we approach the bifurcation. 

Unfortunately this analysis cannot replace that of the previous section, since we have 
not explicitly found a normalisable negative mode of d22l) for a perturbation about the AdS- 
Taub-Bohr instanton. It is possible to find the linearised perturbation of the metric at this 
point using our analysis, but this perturbation is not normalisable. It is however not in the 
Transverse, Tracefree gauge and it is possible that in the appropriate gauge this perturbation 
does give a normalisable mode. 

5. The Continuous Spectrum of A L 

We have been so far concerned with looking for eigenvalues of A L on the AdS-Taub-Nut 
and -Bolt instantons. In the case of an elliptic operator on a non-compact manifold, the 
spectrum may also include a continuum of approximate eigenvalues. As an example, the time 
independent Schrodinger operator for the Hydrogen atom is a Laplacian type operator acting 
on R 3 / {0}. As is well known, this has a countable set of negative eigenvalues corresponding to 
the bound states. The spectrum however includes a continuum of positive "eigenvalues" which 
correspond to free particles moving in the Coulomb potential. These are not true eigenvalues 
since they do not correspond to normalisable eigenfunctions (wavefunctions). The sense in 
which they may be considered part of the spectrum is given below. We shall show that the 
continuous spectrum of Al acting on any of the three spaces considered above includes the 
ray {A e R : A > 9/4Z 2 }. 

5.1. A Little Functional Analysis 

We can think of the space of gauge fixed, finite metric perturbations as a Hilbert space, defined 




(37) 



by: 



H = {h ab E T*M ® s T*M : V a h ab = 0, h 



a 



a 







h ab h ab dV < oo}, 



(38) 



M 



with L 2 



inner product given by 




(39) 
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A L is then a linear operator from 7Y 2 to Hq, where by Hi we mean 

tii = {h ab G H : d Cl . . . d Ci h ab exist and are continuous}. (40) 

Hi are dense subsets of H with respect to the norm topology, where the norm is defined as 
usual by: 

\\h\\ 2 = (h,h) = [ h ab h ab dV. (41) 
Jm 

A linear operator A : H —> H is said to be bounded if 3 A; such that 

\\Ax\\ < k\\x\\ , VxeH, (42) 

an operator which is not bounded is unbounded. We define the spectrum of A, a (A) as 
follows: 



A ^ o-(A) -<=>- the following properties hold: 

(i) (A - A/)" 1 exists, 

(ii) (A - A/)" 1 is bounded, 

(iii) (A — A/) -1 is densely defined (i.e. defined on a dense subset of H). 



The set of A such that (A — A/) -1 does not exist is called the point spectrum of A, Pa (A), the 
set of A such that (A — A/) -1 ) is unbounded is called the continuous spectrum of A, Ca(A). 
We shall not discuss the third condition. 

Lemma 5.1. Given A ^ Pa(A), if there exists a sequence x n G H such that 
\\(A-\I)x n \\ 

ii — 1| ► 0, as n — >• oo, 

II "^n || 

then A G Ca(A). 

Proof. Define y n = (A — \I)x n , by assumption as n — > oo 

IkH n . \\{A-\I)^y n \\ 

U =^ — ► oo, 



\\(A-\I)-iy n \\ \\y n \ 
but this clearly cannot occur if (A — XI)^ 1 is bounded, since if this is the case, 3 k such that 
||(A-AI)-V|| ^ k 

WVnW 

Thus (A — A/) -1 must be unbounded, hence A G Ca{\). □ 
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We now consider the continuous spectrum of A L . We will restrict attention for the moment to 
perturbations of the form considered in section 3. These are generated by the h 00 component, 
which has to satisfy equation d25l) . We first define the functions p n (x) by: 



1 



i 



(x— n) 2 



(x— n— l) 2 







if r < n 
dx ifn<r<n + l 
if r > n + 1 



(43) 



with K chosen so that p n {n + 1) = 0. These functions are everywhere smooth on the real 
line. For A > 9/4Z 2 we define f(r) to be the solution of (1251 which gives a finite contribution 
to (h, h) at r+. This generates a perturbation which satisfies 

(A L - XI)h ab = 0. 

However this h ab is not an element of H since as r — >• oo we have 



f(r)~pr 11//2 cos(a logr) + qr 11/2 sin(a logr 



11/2 



(44) 



where a = ^V4/ 2 A — 9 is a constant. We find that for large r, counting powers of r gives: 



r 

so that (/i, fa) does not converge. We define f n (r) by 
f n (r) = p n (logr)/(r). 



(45) 



(46) 



Setting /i o = /n generates a perturbation which we call h n ah via the constraint equations and 
we claim that 



(A L - XI)h 



ab 



\h n I 

1 11 ab I 



as n 



oo. 



It is easily seen that \\h 



n 
ab 



oo since the integral in ||/i" 6 || is bounded below by the integral 



for \\h ab \\ cut off at r = e n . It suffices then to show that 



(A L - \I)h n ab 



is bounded as 



n — >• oo. Clearly (A L — \I)h™ b = on {r + < r < e n } U {e n+1 < r} so we need to estimate 



XI )K 



ab 



A, 



A/ fa n 



(i6 



(A L - A/)/i r 



(47) 



for large n. In this limit, we can use the large n asymptotic expansions to estimate the leading 
order behaviour of this term. Substituting these in, we find that: 



(A L - XI)h n 



K 



(Lf n ) 2 r w dr. 



with 



Lf n = r 2 f: + l2rf n + (2% + l 2 X)f n 



(48) 



(49) 
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the leading order expansion of (1251) and where /„ is defined by: 

fn( r ) = (pr~ u / 2 cos(a logr) + gr~ n / 2 sin(a log r))p n (logr). (50) 
Under a change of variables x = logr the integral in d48b becomes: 

pe n+1 pn+1 

/ (Lf n ) 2 r w dr = / {p'^(x)(psmax + qcosax) + 2ap' n (x)(qsmax — pcosax)} dx 

J e n Jn 







{Po(y)(P s ^ na (y + n ) + <?cosa(y + n)) 

+2apQ(|/)(gsin + n) — pcosa(?/ + n))} 2 dy, (51) 
since p n (x) = po(x — n). Now for all n the integrand is bounded since 

|Po(*)l <6, 
IPo^l <60, 
thus there exist constants C and n such that 



(A L - XI)h r 



<C, V n>h. (52) 



We can now apply lemma l5~T1 to the operator A L and we find that for A > 9/4/ 2 , (A L — A/) -1 
is unbounded, hence {A > 9/4Z 2 } C Ca(A L ). 

This result makes no use of the value of m or the location of the zeroes of A{r), so it holds 
for all three of the instantons considered here. It is also possible to perform this analysis for 
the other SU(2) invariant modes of the perturbation and we find exactly the same condition 
on A in order that it be in the continuous spectrum. 

6. Conclusions 

We have studied the semi-classical stability of the AdS-Taub-Nut instantons and the two 
branch family of AdS-Taub-Bolt instantons and we have found a negative mode of (|2"2"1) for 
the AdS-Taub-Bolt~ instanton, implying that this instanton is unstable. For the AdS-Taub- 
Nut and AdS-Taub-Bolt + instanton we have argued that they are (linearly at least) stable. We 
have also justified this by considering a family of regular metrics, not necessarily satisfying 
Einstein's equations, which contains these instantons as special cases. This gives an intuitive 
sense of how the saddle-node bifurcation in the S — p plane arises, as well as the Hawking- 
Page type bifurcation which occurs. We have also found that for all three instantons the 
continuous spectrum includes the ray {A > 9/4Z 2 }. 
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The following are a set of differential equations governing the J = M = K = piece of the 
tensor harmonic decomposition of the equation 

— V e V e /i a fe — 2R ac bdh cd = Xh a b, 

for the metric 

ds 2 = -j^dr 2 + 4N 2 A{r)a 2 3 + 5(r)(a 2 + a 2 ). 
Ayr) 

We suppress the r dependence of A(r) and B(r) and use a prime to denote differentiation 
with respect to r. These are taken from lfT2l . with some typos corrected. 
11 Component 

a{r)h'[ x + b(r)h' n + c(r)h 00 + d(r)h u + e(r)h 22 + J(r)h 33 = -Xh u , (A.l) 



with 



a(r) = A, 

b(r) =A'- A—, 

AA2r>l2 

c(r) = - 2A(A'B') + IfL^- _ 4A 2 E", 

B 

2 A'-B' 1 . f B'\ 2 AN 2 A 1 



1 4iV 2 A I (B' 

J{r> 8N 2 A B' 
22 Component 

a ( r )h 2 2 + K r ) h 22 + c(r)h 00 + d(r)h u + e(r)h 22 + f(r)h 33 = -Xh 2 2, (A.2) 

with 

a(r) = a(r), d(r) = e(r), 

b(r) = b(r), and e(r) = d(r), 

c(r) = c(r), /(r) = J(r). 

Where the overbar denotes the coefficients for (IA.1I) . 
00 Component 

a(r)/& + &(r)/i' 00 + c(r)/i 00 + d(r)X + f(r)h 33 = -Xh 00 , (A3) 

with 

a(r) = A, 

4 R' 

6(r) =3A' + ^-, 
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, . A' 2 A'B' . (B'V A „ 

A'B' B' 2 B" 

d(r) = — 



8AB 2 AB 3 4.B 2 ' 

f( , J? A»_ 

J{r> 32N 2 A 3 16N 2 A 2 ' 
and X = hn + h 22 . 
33 Component 

a(r)/& + b(r)h' 33 + c(r)h 00 + d{r)X + f(r)h 33 = -Xh 33 , (A.4) 

with 

a(r) = A, 

6(r) = - A' + — , 

c(r) = 8iV 2 A4' 2 - 16iV 2 A 2 A", 

. 16iV 4 A 2 2AT 2 AA' J B' 
d (r) = — W3 W —, 

t „ A' 2 A'B' AN 2 A 
/W= - A + 2A~^-^ 
andX = hn + h 2 2- 

Constraint Equations 

The above system of equations is consistent with the transverse and tracefree conditions 
on h ab . This gives us the constraint equation 

Ah 00 + -L{h n + h 22 ) + j^jh, = 0, (A.5) 

coming from the traceless condition h a a = 0, and 

(3 AB'\ B' A' 

2 A ' + ~B~) h °° ~ SB 2 ' {ku + k22) ~ ^PA 2 ^ = °' (A - 6) 



which comes from the transverse condition V ' a h 



ab 
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